Effective Abelian and non-Abelian gauge potentials in cavity QED 
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Cavity QED models are analyzed in terms of field quadrature operators. We demonstrate that in 
such representation, the problem can be formulated in terms of effective gauge potentials. In this 
respect, it presents a completely new system in which gauge fields arise, possessing the advantages of 
purity, high control of system parameters as well as preparation and detection methods. Studying 
three well known models, it is shown that either Abelian or non-Abelian gauge potentials can 
be constructed. The non-Abelian characteristics are evidenced via numerical simulations utilizing 
experimental parameters. 
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Introduction. - Gauge fields naturally arise when de- 
scribing subatomic interactions. Ranging from classical 
electromagnetism to quantum Hall systems ^ and 
more recently cold atoms in optical lattices [3j], gauge 
theories have deepen our understanding for fundamental 
processes in AMO physics. The simplest and most famil- 
iar example is found when considering a charge particle 
in an electromagnetic field. In this case, the gauge the- 
ory is Abelian, as the vector components of the gauge 
field mutually commute. For non-Abelian gauge fields, 
on the other hand, the vector field components are non 
commuting operators. In general, time-ordering then be- 
comes important, leading to novel phenomena. With the 
experimental progress in especially AMO physics, non- 
Abelian gauge structures have drawn great interest in 
recent years. 

Wilczek and Zee showed that by adiabatically chang- 
ing a Hamiltonian, possessing degenerate states, effective 
non-Abelian gauge potentials can be achieved In- 
deed, the appearance of gauge potentials for adiabatic 
evolution had been demonstrated some years prior by 
Mead and Truhlar studying molecular systems [5|. Here, 
adiabaticity is a result from separation of fast electronic 
and slow nuclear motions in the molecule, i.e. the Born- 
Oppenheimer approximation. Thus, the gauge potential 
derives from intrinsic spatial evolution [y] and not from 
explicit time-dependence as utilized in [J]. Similar situ- 
ations emerge for cold atoms interacting with spatially 
varying light fields 0] . The center of mass motion of the 
atom induces an effective gauge potential, and by consid- 
ering a coupled four-level tripod atomic system, Ruseckas 
et al. presented a model which exhibits a non-Abelian 
gauge structure ^. The advantage of this cold atom 
model is the high controllability of system parameters as 
well as efficient preparation and detection methods. 

In this Letter we present a completely different sys- 
tem in which effective gauge potentials appear, with the 
same assets as for the cold atom model, i.e. purity, co- 
herent control of system parameters, preparation, and 
detection. In particular, we consider a single few-level 



atom (could be a true atom or a quantum dot repre- 
senting an artificial atom) interacting with one or two 
quantized cavity modes. By expressing the field in terms 
of its quadrature operators, the structure of the Hamilto- 
nian becomes similar to the ones encountered in molecu- 
lar models. In other words, an adiabatic diagonalization 
renders effective gauge fields. This is demonstrated by 
analyzing three well known models, where two are en- 
dowed with non-Abelian properties. 

Since two decades, cavity quantum electrodynamics 
(QED) with single or few atoms has delivered some of the 
most striking experimental results on pure quantum phe- 
nomena 9]. Among others, entanglement generation [lO|, 
the quantum measurement problem and the quantum- 
classical transition 1J| , and verification of the graininess 
of the quantized electromagnetic field Cavity QED 
has attracted even more interest in recent years due to 
the realization of coherent coupling of single quantum 
dots [3] or Bose-Einstein condensates [31 to a cavity 
mode. These experiments pave the way for the possibil- 
ity of reaching a super-strong coupling regime of cavity 
QED. 

The general form of our Hamiltonian reads 



H = Hf + Ha + Hi 



(1) 



where 



N 



H, = Y,E,\j){j\, 



(2) 



Here, ujk is the /c'th field mode frequency, aj. (ak) the cre- 
ation (annihilation) photon operator of mode fc, Ej the 
energy of atomic level j, N the number of atomic states, 
d^^^f the dipole operator for atomic transition j, and ^(x) 
is the electric cavity field. In the dipole approximation 
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we set X = 0, and the field can be written as 



E 



ak - a\ 



(3) 



where £k is the polarization vector for mode k and £k 
the corresponding field amplitude. Moreover, the compo- 
nents of the dipole moment become d!£^ — ~e\j){j\a\j + 
1)0 + 1| + h.c. with a = x,y,z, e the electron charge, 
and h.c. is the Hermitian conjugate. For our purpose, 
it is convenient to express the Hamiltonian in its field 
quadrature operators 



Pk = 



ak - al 



(4) 



obeying the canonical commutation relations; [Xi,,Pfe/] — 
iSkk'- The field quadrature operators ^ are easily 
measured experimentally. Indeed, the ENS group of S. 
Haroche recently presented experimental results were the 
full phase space distribution of a cavity mode was as- 
sessed 15|. In terms of (|4]), the Hamiltonian ([1]) takes 
the form 



H = 



N 



k 
N 



(5) 
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where gjk is the effective atom-field coupling between the 
transition j and the mode k. The general form of the 
Hamiltonian ([5]) works as a starting point to analyze more 
specific cavity QED models. 

Rabi model. ~ The simplest non-trivial situation con- 
siders a single cavity mode interacting with one atomic 
transition. Experimentally, such idealized situations are 
accessible by utilizing high-Q cavities with well resolved 
resonance frequencies Wfc, see for example Refs. 0, [13, 



111 . \12\. For this model we leave out the indices j and 
k. Introducing the Pauli matrices (Jx = |2)(1| -I- |1)(2|, 
ay = -^|2>(1| +^|1)(2|, and a, = |2)(2| - |1)(1|, and 
without loss of generality assuming that d ■ e is purely 
real, we obtain the Rabi Hamiltonian 

+ — I + -Y^^ ~ f^g^xP- (6) 

Here, g is the vacuum Rabi frequency, giving the effective 
atom-field coupling, and we have chosen our zero energy 
such that El =~Ml/2 and E2 = Ml/2. 

We note that the Hamiltonian can be rewritten as 



HRab 




(7) 



where we have introduced the scaled gauge potential A = 

2 

^(Jx and the scalar potential $ = — ftfj. Under a unitary 



transformation of the state vector 

■^{x,t)^u\x,t)-^{x,t), 



(8) 



the gauge and scalar potential transforms accordingly 
A U\X,t)AU{X,t) - u\X,t)-^U(X,t), 

(9) 

<! ^ u\x, t)m{x, t) - ihu\x, t)—u{x, t). 

The above equation demonstrates the gauge structure of 
the effective potentials. 

It should be observed that the appearance of the ef- 
fective gauge potentials results from the non-stationary 
dynamics of the quantized cavity field, and not from adi- 
abatic particle motion as in earlier works on cold atoms 
and molecular physics [^, 0, S Hi- Furthermore, note 
that by applying the rotating wave approximation, Hjiabi 
turns into the solvable Jaynes-Cummings Hamiltonian 
[3] that has served as workhorse in the field of quan- 
tum optics ever since it was introduced [17|. Imposing 



such an approximation would make the present analysis 
less transparent since the atom field coupling would then 
depend on both P and X rendering complicated gauge 
potentials. In addition, for current quantum dot cavity 
QED systems, application of the rotating wave approxi- 
mation is not always justified. 




Figure 1: Effective APSs of the bimodal Rabi Hamiltonian 
(|10p in the g > y/ujQ case, where the lower APS attains the 
sombrero shape. The conical intersection is located at the 
origin. Pi = P2 = 0. 

Bimodal Rabi model. - In order to acquire non-Abelian 
gauge potentials, additional degrees of freedom must be 
included. This is easily accomplished in our cavity QED 
setting by considering a bimodal cavity field [3]. Hence, 
the atom interacts simultaneously with two cavity modes. 
It is worth noting that experiments on bimodal cavities 
have been successfully demonstrated [31 . The simplest 
extension of the Rabi model is to keep the two-level struc- 
ture of the atom and simply add one additional cavity 
mode. Recently, it was shown that the corresponding 
cavity QED model has a Jahn- Teller structure 20]. For 
appropriate choices of polarizations and atomic dipole 



3 



moments, one obtains the E x e Jahn- Teller Hamiltonian 
which has been thoroughly studied in especially molec- 
ular physics [2l|. The model possesses a conical inter- 
section in which the adiabatic potential surfaces become 
degenerate. Encircling the intersection brings about a 
non-zero geometrical phase characterized by a gauge po- 
tential. We note that geometrical phases have been dis- 
cussed in terms of cavity QED 2^ 22 1. These, however, 
do either consider the Abelian situation or Berry phases 
originating from a time-dependent Hamiltonian and con- 
sequently the gauge field derives from external driving 
rather than from intrinsic dynamical evolution. 

The E X e cavity QED Hamiltonian is given by f20| 



0/{a) 



BR 




(10) 



Here, we have assumed equal mode frequencies and equal 
atom- field strengths. Contrary to Jahn- Teller models 
encountered in molecular physics, the conical intersec- 
tion appears in momentum space rather than in po- 
sition. In the condensed matter community, a cou- 
pling as the one in Eq. (|10p is usually said to be on 
Rashba form |23|]. The adiabatic potential surfaces 
(APS), defined as V;=^(Pi,P2) - huj {P^ + P^) /2 ± 
+ {Pi + P^), are envisaged in Fig. [1] Due 
to the non-zero Q, the conical section becomes avoided. 
The lower APS has a sombrero shape whenever g > \/Quj, 
otherwise a global minimum at the origin is exhibited. It 
follows directly from the form of the Hamiltonian pO| . 
that the vector and scalar potentials read 



id 
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(11) 



and since 
Abelian. 

In a recent work 
Einstein condensate 



[^1,^2] ^ 0, the gauge potential is non- 



2J], considering a spinor Bose- 
a proposal for detecting non- 
Abelian characteristics was put forwarded. The idea is to 
initialize a state, such that it is given by a fairly localized 
wave packet with non-zero average momentum, and then 
boost the wave packet either clockwise or anti-clockwise. 
Due to the non- Abelian structure, the two paths will ren- 
der different dynamics despite the polar symmetry of the 
problem. To demonstrate this, we chose the atom to be 
initially in a superposition state (|2) — |1)) /\/2 and the 
two fields to be in coherent states 



1 



.^-{'^■P.oy^-^X,oP.^-iP.-P.or/2^ i = 1, 2, 



(12) 

where Pio and Xio are initial average momentum and 
position respectively. These are related to the initial co- 
herent state amplitude aio = (X^o -I- iPiQ)/\/2. For the 



S -0.2 




Figure 2: Time evolution of the atomic inversion for clockwise 
(a) and anti-clockwise (b) propagation. Solid lines display the 
ideal case of no losses, while for the dotted lines both cavity 
and atomic losses are included. The non-Abelian structure 
of the model is clearly visible. The system parameters are 
Ao = 2, P20 = 0, Xio = 0, X20 = ±5, Q./2-K = 6.9 GHz, 
uj/2ti = 5.7 GHz, g/2-n: = 105 MHz, 'y/2n = 1.9 MHz, and 
k/2tt = 250 kHz. 



situation at hand, we set Pio > 0, P20 — 0, Xiq = 0, 
and X20 7^ 0. In momentum representation, this gives 
an initial wave packet mainly at the lower APS located 
at the positive Pi-axis and with an initial velocity per- 
pendicular to this axis. For X20 > 0, the wave packet 
sets off clockwise around the origin, while A"2o < re- 
sults in anti-clockwise evolution. Our numerical simula- 
tion utilizes the split-operator method, which gives the 
time-evolved wave packet at any instant of time. Solid 
lines of Fig. [5] presents the results for the atomic inver- 
sion, W{t) — P2{t) — pi{t) where Pi{t) is the probability 
of finding the atom in the state \i) at time t, for clockwise 
propagation (a) and for anti-clockwise evolution (b). The 
final time corresponds to approximately six roundtrips 
around the conical intersection. We point out that the 
atomic inversion can be experimentally measured up to 
a few percent accuracy ^ . The discrepancy between the 
lines of (a) and (b) , evident even at very short time scales 
(~ns), is a manifestation of the underlying non-Abelian 
character. 

A deeper understanding is obtained by studying the 
evolution of the fields in phase space. As for a classical 
harmonic oscillator, both fields will encircle the origin of 
phase space with, for the example of Fig. [21 radii approx- 
imately 2 and 5 respectively. The effective magnetic field 
deriving from the gauge potential renders a momentum 
dependent force. The momentum is in general differ- 
ent for clockwise and anti-clockwise evolution, and hence 
this effective Lorentz force acting on the field distribu- 
tions implies slightly deviating trajectories in phase space 
for the two cases. This difference is the origin for the 
momentum-dependence of the atomic inversion seen in 
Fig. [2 

In realistic situations, both cavity and atomic losses 
come into play. We schematically take these losses 
into account by consider time evolution of the ef- 
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fective non-hermition Hamiltonian H, 
IK ^a|ai + 0,202^ — i7|2)(2|. Here k and 7 are the photon 
decay rate and the atomic spontaneous emission rate re- 
spectively. The results of such numerical simulation are 
given by the dotted lines of Fig. [21 The parameters used 
in the examples of Fig. [2] have been chosen according to 
the quantum dot cavity QED experiments of Ref. 
Due to the short interaction times of Fig. [2l losses play a 
minor role. The importance of losses becomes apparent 
only for longer time scales. 

Bimodal A model. - Our final example considers a 
A-atom with two lower meta stable states |1) and |2) 
coupled to an excited state |3) via two cavity fields 1 
and 2 respectively. Numerous theoretical works on this 
system have been put forward , and Ref. [l^l presents 
some experimental studies of A cavity experiments. 

The Hamiltonian for this system reads 




(13) 



where we have assumed that the effective atom-field cou- 
plings are real and the same between the two transitions. 
For degenerate lower atomic states Ei = E2, the three 
APSs possess a Renner- Teller intersection [26] between 
the lower and the middle surface and another one be- 
tween the middle and the upper surface. These intersec- 
tions are characterized by that the tangents of the two 
surfaces are the same at the degenerate point. Moreover, 
the Berry phase acquired by encircling a Renner- Teller 
intersection in position space vanishes. Here, however, 
the intersection is in the momentum space giving rise 
to a non-Abelian Berry phase in position space, where 
we have checked that for Ei = E2 the diagonal terms 
of the corresponding 3x3 geometric phase matrix are 
zero but some of the off-diagonal terms are indeed non- 
zero. If the two lower atomic states are not degenerate 
El ^ E2, both Renner- Teller intersections split into two 
non-avoided conical intersection. 

As in the previous bimodal example, the gauge poten- 
tial is non-Abelian; 



(^l,A2) = ^(A4,A6) 
UJ 




(14) 
|3)(1| + 

+ |2)(3|, and Ag = (|1)(1| + |2)(2| - 



with Xj being the Gell-Mann matrices (A4 



|1>(3|, Ae = |3)(2| 
2|3)(3|)/V3). 

Concluding remarks. - We have demonstrated how 
cavity QED models provide novel systems exhibiting ar- 
tificial gauge potentials. This derives from the quantized 
motion of the cavity fields; the dynamics of the fields in 
phase space. Moreover, using numerical simulations we 



showed that non-Abelian characteristics should be de- 
tectable under realistic conditions. From this rather un- 
usual approach to cavity QED, presented in this Letter 
and earlier in [3, H^l j it is clear that these models are 
very rich. We are at the moment analyzing the prospects 
of achieving zitterbewegung or Hall effects with cavity 
QED setups [27[ . We hope that the current contribution 
will encourage experiments along these lines. 

We wish to thank Prof. Erik Sjoqvist for fruitful dis- 
cussions, and JL acknowledge support from the MEG 
program (FIS2005-04627). 
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